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Abstract 

The effects of both kinematical and dynamical SU(6) breaking on the nucleon elastic 
form factors, G E (Q 2 ) and G^(Q 2 ), are investigated within the constituent quark model 
formulated on the light-front. The investigation is focused on G E (Q 2 ) and the ratio 
G P m(Q 2 )/Gm(Q 2 ), which within the SU(6) symmetry are given by G E (Q 2 ) = and 
(*m(Q 2 )/ ^%i{Q 2 ) = — 3/2, respectively. It is shown that the kinematical SU(6) 
breaking caused by the Melosh rotations of the quark spins as well as the dynamical 
817(6) breaking due to the mixed-symmetry component generated in the nucleon wave 
function by the spin-dependent terms of the quark-quark interaction, can affect both 
G n E {Q 2 ) and G p M {Q 2 )/G n M {Q 2 ). The calculated G E (Q 2 ) is found to be qualitatively 
consistent with existing data, though only ~ 65% of the experimental neutron charge 
radius can be explained without invoking effects from possible non-vanishing sizes of 
the constituent quarks and/or many-body currents. At the same time the predictions 
for the magnetic ratio G P M (Q 2 ) / G^^Q 2 ) turn out to be inconsistent with the data. 
It is however shown that the calculations of G^j(Q 2 ) based on different components 
of the one-body electromagnetic current lead to quite different results. In particular, 
the calculations based on the plus component are found to be plagued by spurious 
effects related to the loss of the rotational covariance in the light-front formalism. 
These unwanted effects can be avoided by considering the transverse y-component of 
the current. In this way our light-front predictions are found to be consistent with the 
data on both G E (Q 2 ) and G p M (Q 2 )/G r l I (Q 2 ). Finally, it is shown that a suppression of 
the ratio G P E (Q 2 ) / G P M (Q 2 ) with respect to the dipole-fit prediction can be expected in 
the constituent quark model provided the relativistic effects generated by the Melosh 
rotations of the constituent spins are taken into account. 



1 Introduction 



The nucleoli elastic electromagnetic (e.m.) form factors contain important pieces of infor- 
mation on the internal structure of the nucleon, and therefore an extensive program for their 
experimental determination is currently undergoing and planned at several facilities around 
the world Ji]]. In this work we focus on the nucleon Sachs form factors, defined as |§ 

G N E (Q 2 ) = f^q 2 ) - t^W) . 

g n ,aq 2 ) = F?m+F 2 N m (i) 

where Ff(Q 2 ) [FJ?(Q 2 )} is the Dirac [Pauli] form factor appearing in the covariant decom- 
position of the (on-shell) nucleon e.m. current matrix elements, viz. 

< N(p', s')\j» m (0)\N(p, s) >= u(p', s') {F^Q 2 )^ + F?{Q 2 ){ia^q v /2M)}u{ P , s) (2) 

with Q 2 = —q ■ q and M being the squared four-momentum transfer and the nucleon mass, 
respectively. As it is well known (cf., e.g., Ref. [[J), the spin-flavour SU(6) symmetry 
predicts G n E (Q 2 ) = and G p M (Q 2 )/G n M (Q 2 ) = -3/2. 

The term (-Q 2 /4M 2 )F 2 W (Q 2 ), appearing in the definition of G% (Q 2 ) (Eq. (0)), is usually 
referred to as the Foldy contribution and it is of relativistic origin. Its contribution to the 
neutron charge radius is known to yield almost totally the experimental value of the latter 
(( r )n,ex P = —0.113 ± 0.005 fm 2 from Ref. Q). This result has been viewed || as an 
indication of the smallness of the intrinsic charge radius related to the neutron rest-frame 
charge distribution. Very recently || however the interpretation of the neutron charge radius 
as arising from its internal charge distribution has been addressed again, arguing that, going 
beyond the non-relativistic limit when the Foldy term firstly appears, the Dirac neutron 
form factor F™(Q 2 ) receives a relativistic correction that cancels exactly against the Foldy 
term in G%(Q 2 ). Such a statement was inferred from the observation that the well-known 
phenomenon of zitterbewegung, which produces the Foldy term, cannot contribute to the 
neutron charge radius, because the latter has zero total charge JJ. 

The result of Ref. || has been recently derived from a relativistic calculation of the 
nucleon elastic e.m. form factors, performed within the light-front (LF) formalism and 
the constituent quark (CQ) model 0. Namely it has been shown that the zitterbewegung 
approximation of Ref. || corresponds to neglect the initial transverse motion in the Melosh 
rotations of the CQ spins. In Ref. the predictions of both the non-relativistic (NR) 
limit and the zitterbewegung (ZB) approximation have been carried out assuming a pure 
S'f/(6)-symmetric nucleon wave function. In both cases it was found that G^(Q 2 ) = and 
G\i(Q 2 ) / G n M {Q 2 ) = —3/2. While the latter is consistent with experimental data, the former 
prediction is well known to be at variance with existing data. Therefore, both in the NR 
limit and in the ZB approximation the only way to produce a SU(6) breaking can be of 
dynamical origin, i.e. related to the presence of a mixed-symmetry component generated in 
the nucleon wave function by the spin-dependent terms of the effective quark-quark potential. 
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We point out that spin-dependent forces are important at short interquark distances, and 
their strengths are almost fixed by the spin splittings in the hadron mass spectroscopy (cf., 
e.g., Refs. and §). 

However, in Ref. || it has been shown that: 

• the effects of the CQ initial transverse motion in the Melosh rotations can be neglected 
to a good approximation only when the average value of the transverse momenta, (p±), 
is much smaller than the CQ mass m. However, in QCD both m and (p±) are expected 
to be at least of the order of the QCD scale, Aq C d ~ 300 MeV. Moreover, in quark 
potential models, like those of Refs. 0] and ||, (p±) turns out to be much larger than 
m, because of the high momentum components generated in the hadron wave functions 
by the short-range part of the effective CQ interaction (cf. Ref. [[|); 

• the Melosh rotations produce a kinematical SU(6) breaking. Indeed, such rotations, 
being momentum and spin dependent, produce a re-coupling of the CQ orbital angular 
momentum and spin; in other words, the nucleon LF wave function cannot be any more 
expressed as a product of a spatial part times a spin-isospin function, and therefore it 
cannot be SU(6) symmetric. The inclusion of the effects of the CQ initial transverse 
motion lead unavoidably to G%{Q 2 ) ^ and G P M {Q 2 ) / G'l^Q 2 ) ^ —3/2, even when 
the dynamical SU(Q) breaking, generated by the spin-dependent components of the 
effective CQ interaction, is neglected. Moreover, it has been shown that LF relativistic 
effects may contribute significantly to G%{Q 2 ) when (p±) > m, and in particular the 
calculated neutron charge radius can reach ~ 40% of its experimental value; 

• the SU(6) breaking associated to the Melosh rotations heavily affects also the ratio 
G P M {Q 2 ) / G^iQ 2 ) , leading however to a significative underestimation of the experimen- 
tal data. 

The aim of the present work is twofold: i) to include the effects of the mixed-symmetry 
wave function both in the NR limit and in the ZB approximation as well as in the full LF 
approach, extending in this way the analysis of Ref. ||; ii) to show that LF calculations of 
Gm(Q 2 ), based on different components of the one-body e.m. current, lead to quite different 
results, because of spurious, unphysical effects related to the loss of the rotational covariance 
in the light-front formalism, while the same does not occur in case of the nucleon charge 
form factor G%(Q 2 ). 

It will be shown that: i) the predictions of the NR limit and the ZB approximation 
underestimate significantly the data on G^(Q 2 ), even when the effects of the mixed-symmetry 
wave function are considered, and ii) once spurious effects are properly avoided, the LF 
predictions are consistent with the data on both G%{Q 2 ) and G P M (Q 2 ) / G^iQ 2 ) , though only 
~ 65% of the experimental neutron charge radius can be explained without invoking effects 
from possible non- vanishing CQ sizes and/or many-body currents. This means that the 
dynamical SU(6) breaking, predicted by quark potential models based on the hadron mass 
spectroscopy, appears to be consistent with the SU(6) breaking observed in the nucleon 
elastic e.m. form factors, provided relativistic effects are properly taken into account. 
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Finally, we will address also the issue of the interpretation of the recently observed [HJ de- 
viation of the ratio fj, p G p E (Q 2 ) /G P M (Q 2 ) from the dipole-fit expectation n p G p E {Q 2 ) / G P M {Q 2 ) = 
1, where /i p = G P M (0) is the proton magnetic moment. It will be shown that a suppression 
of the above mentioned ratio from unity can be expected in the CQ model provided the 
relativistic effects generated by the Melosh rotations of the CQ spins are taken into account. 
The results of the calculations obtained using the nucleon wave functions arising from two of 
the most sophisticated quark potential models 0, § , are found to compare quite favourably 
against the recent data [10] from Jefferson Lab [J LAB). 



2 The nucleon light-front wave function 

Let us briefly recall the basic notations and the relevant structure of the nucleon wave 
function in the LF formalism. As it is well known (cf., e.g., Refs. [O, [T3, [TBI]), the nucleon 
LF wave function is eigenstate of the non-interacting (LF) angular momentum operators j 2 
and j n , where the unit vector n = (0, 0, 1) defines the spin quantization axis. The squared 
free-mass operator is given by 



Ml 



£(l^l 2 + ™ 2 )/& 



(3) 



i=i 



where 



t _ Pi 
& ~ p+ 



Pi± - &P± 



(4) 



are the intrinsic LF variables. The subscript _L indicates the projection perpendicular to 
the spin quantization axis and the plus component of a 4- vector p = (p°,p) is given by 
p+ = pO _|_ ^ . p*. g na Hy p = (p+ ; p ± ) = p 1 4- p 2 4- p 3 is the nucleon LF momentum and pi 
the quark one. In terms of the longitudinal momentum k in , related to the variable £j by 



h- — 

ru tn — 



&M - (\k t± \ 2 + m 2 ) /&M 



the free mass operator acquires a familiar form, viz. 

3 3 



(5) 



(6) 



i=l 



i=l 



with the three-vectors ki defined asQ 

ki = (ki±, ki n ) 

"Note that ki are internal variables satisfying k\ + ki + £3 = 0. 



(7) 



4 



Disregarding the colour variables, the nucleon LF wave function reads as 



E1E2E3 



ut(t E ({^I^IMXfe^lx?) (8) 

V JWo4l4243 {n} 

where the curly braces { } mean a list of indexes corresponding to % — 1, 2, 3, and (r^) is 
the third component of the CQ spin (isospin). The rotation 7U, appearing in Eq. (Q), is the 
product of individual (generalized) Melosh rotations, viz. 



7^ = n^(w. 



where 



m) 



(9) 



(10) 



with a being the ordinary Pauli spin matrices. 

In what follows we neglect the very small P- and D- waves of the nucleon (see later on), 
and therefore we limit ourselves to the following canonical (or equal-time) wave function 
(corresponding to a total orbital angular momentum equal to L — 0): 



{{k;viTi}\xN N ) 



C ({^})-< TJV ({^» 
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vntn 1 



+ 
+ 
+ 



w A (k,p) 



V2 



01 

V N T N 



{{nn}) - ({^}) 



where Ws{k,p), Ws' s (k,p), ws> a (k,p) and WA{k,p) are the completely symmetric (5), the two 
mixed-symmetry (S") and the completely antisymmetric (A) wave functions, respectively. 
The variables /c and p are the Jacobian internal coordinates, defined as 



k 



h - k 2 



P 



2k 3 - (hi + k 2 



(12) 



with ki given by Eq. (0). Finally, the spin-isospin function $^ 2 ^ 2 ({i/jTj}), corresponding to 
a total spin (1/2) and total isospin (1/2), is defined as 

$5S 2 ({^}) = E(\^\S 12 M S ) (S 12 M s ~v 3 \±v N ) ■ 

M s 
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X](^ r l^ r 2l T 12^T) (Ti 2 M T ~T 3 \-T N ) 



(13) 



where Su (T12) is the total spin (isospin) of the quark pair (1,2). The normalizations of the 
various partial waves in Eq. (ITT]) are 



dk dp ws{k,p) 



dk dp ws> (k,p) 



dk dp Wji(k,p) 



Ps 



dk dp ws' a {k,p) = Ps'/2 



Pa 



with P s + P s > + P A = 1. 

Generally speaking, the nucleon LF wave function 
of a mass equation of the form 



tf^) (see Eq. 



(14) 



is eigenfunction 



{M + V}\^ N N ) 



(15) 



where Mq is the free-mass operator @ and V is a Poincare-invariant interaction term. Using 
the shorthand notation \^ v ^) = 7V\x'n') 1 where \xn) ls the nucleon canonical wave function 
(see Eq. ([□])), one gets 



.i=i 



m 2 + 



h\ 2 + v \\x7) = m\x v n n ) 



(16) 



where V = TZV7V is the Melosh-rotated interaction term, while the free-mass Mq is invariant 
under Melosh rotations and it has been expressed directly in terms of the three-vectors ki 
through Eq. (^|). The Poincare-invariance of V implies simply that V has to be independent 
of the total momentum of the system and invariant upon spatial rotations and translations. 
Therefore, one can adopt for V any quark potential model able to reproduce the hadron 
mass spectra, simply interpreting ki as canonical variables and treating the kinetic energy 
operator in its relativistic form. 

In this work we consider one of the most sophisticated quark potential models, based 
on the assumption of the valence + gluon dominance (the OGE model of Ref. 0). The 
baryon (canonical) wave functions are calculated by solving the corresponding three-quark 
Hamiltonian problem through a variational technique, based on the expansion of the wave 
function into the harmonic oscillator basis (see appendix A for more details). The coefficients 
of the expansion and the eigenvalues of the three-quark Hamiltonian are determined via the 
Raleigh-Ritz variational principle. The presence of spin-orbit and tensor components in 
the OGE model give rise to non-vanishing P- and D-waves in the nucleon, respectively. 
The probabilities of the various partial waves in the nucleon turn out be: Ps = 98.04%, 
P s - = 1.70%, P A < 10" 2 %, P P = 0.05% and P D = 0.21%. It can be seen that due to 
the smallness of the spin-orbit and tensor terms of the OGE model (which is required by 



6 



the smallness of the spin-orbit and tensor splittings in the hadron masses) the P- and D- 
waves in the nucleon are weakly coupled to the dominant symmetric S'-wave. Moreover, 
also the antisymmetric wave WA{k,p) exhibits an extremely weak coupling to the dominant 
symmetric S*-wave. Finally, the probability of the mixed-symmetry S' component is of the 
order of few %, being governed by the strength of the spin-spin component of the OGE 
model, which is almost fixed by the N — A(1232) mass splitting. We point out that the same 
basic features are shared also by the nucleon wave function arising from a recently proposed 
quark potential model, based on the exchange of the pseudoscalar mesons arising from the 
spontaneous breaking of chiral symmetry (the chiral model of Ref. |§). As a matter of fact, 
in case of the chiral model the partial wave probabilities are found to be: Ps = 98.73%, 
P s , = 1.27%, P A < 10- 2 %5 




1 2 3 1 2 3 

p (GeV/c) p (GeV/c) 



Figure 1. Canonical CQ momentum distribution in the nucleon, n(p) (Eq. (|l7|)), versus the 
three-momentum p = |p| (see Eq. (|l2D). (a) Thick and thin lines correspond to the nucleon 
eigenfunctions of the OGE [|J and chiral || quark potential models, respectively. Solid lines are 
the results obtained using the full interaction models, while dotted lines correspond to the inclusion 
of their linear confinement terms only, (b) Contributions of various partial waves (see Eq. (p|)) 
to the CQ momentum distribution n(p) obtained within the OGE interaction model. Dot-dashed, 
dotted and dashed lines correspond to S' a -, S' s - and S- waves, respectively. The solid line is their 
sum. 

In Ref. || the wave functions of the nucleon and of the most prominent electroproduced 
nucleon resonances corresponding to the OGE and chiral quark potential models have been 
compared. Here, we have reported in Fig. 1(a) the (canonical) CQ momentum distribution 

^In case of the version of Ref. H the chiral quark potential model does not contain any spin-orbit and 
tensor terms and therefore the resulting nucleon wave function has no P- and D-waves. 
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n(p) in the nucleoli, calculated within the OGE and chiral models neglecting the small P- 
D- and antisymmetric wa partial waves, namely: 



n(p = \p\] 



dVLp dk 



ws{k,p) + w s > a (k,p) + w S ' a (k,p) 



(17) 



From Fig. 1(a) it can clearly be seen that both quark potential models predict a remark- 
able content of high-momentum components, which are generated by the short-range part 
of the quark-quark interaction (cf. also Ref. ||). In case of the OGE model the interaction 
at short-range is mainly due to the flavour-independent central Coulomb-like and spin-spin 
terms arising from the effective one-gluon-exchange, while in case of the chiral model the 
interaction at short interquark distances is governed by the spin-flavour structure charac- 
terizing the exchanges of pions, rj and partially rf mesons. As a matter of fact, the above 
mentioned high-momentum components are completely absent if one switches off the inter- 
mediate boson exchanges and only the linear confining term (dominant at large interquark 
distances) is considered. In Fig. 1(b) we have reported the individual contributions of the 
S-, S' 8 - and S' a - waves to the CQ momentum distribution n(p), calculated within the OGE 
model. It can be seen that both the S' s and S' a components are negligible at low momenta, 
while they may be important at high momenta, as expected from the short-range nature of 
the spin-spin interaction among quarks. 

Table 1. Values of the average CQ transverse momentum (p±) (f[^), the CQ mass radii (r) u (<n 
( |T9| ) and the charge radii of proton and neutron ( f20| ) predicted by the OGE |7j and the chiral |8[ 
quark potential models. 





OGE model 


chiral model 


conf. only 


w/o S' 


with S' 


w/o S' 


with S' 


(p±> (GeV) 


0.33 


0.56 


0.58 


0.60 


0.61 


(r)u (fm) 


0.49 


0.31 


0.32 


0.31 


0.32 


(r)d (fm) 


0.49 


0.31 


0.31 


0.31 


0.30 


(r) p (fm) 


0.49 


0.31 


0.32 


0.31 


0.32 


(r)l (fm") 


0.00 


0.00 


-0.0083 


0.00 


-0.0090 



Finally, we have collected in Table 1 the predictions for few observables which are of 
interest in this work, namely the (canonical) average CQ transverse momentum (p±) 



{px)=J-J o dpp 4 n(p) 
and the (canonical) CQ mass radii (r) u rd) m the proton 



<r> 



\ z f P i=i z 



(19) 
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where fi is the conjugate variable to fej. In terms of (r) u (d) the squared (canonical) charge 
radii of proton and neutron are given by 

(r)l = l(Mr)l-(r)t) 

(r)l = l((r) 2 d -(r) 2 u ) (20) 

From Table 1 it can be seen that: i) both models give rise to a proton with a small size 
(~ 0.3 fm instead of ~ 0.8 fm), and ii) the predicted neutron charge radius has the correct 
sign, but its absolute value is about one order of magnitude less than the experimental value, 
« e , P = -0-113 ± 0.005 /m 2 @. 



3 Calculations of the nucleon elastic form factors 



As in Ref. (cf. also Ref. || 
operator at the CQ level, viz. 



we consider the one-body component of the e.m. current 



3=1 



3 

E 

.7=1 L 



e 3 Yn(Q 2 ) + 



2m 



IKj 



tt(Q 2 



(21) 



where a Ufl = Ifr", 7**], &j is the charge of the j-th quark, k,j the corresponding anomalous 
magnetic moment and f^ 2 )(Q 2 ) ^ s Dirac (Pauli) form factor (with f^ 2 )(Q 2 = 0) = 1). In 
the LF formalism the form factors for a conserved current can be obtained using only the 
matrix elements of the plus component of the current operator and, moreover, for Q 2 > 
the choice of a frame where q + = allows to suppress the contribution of the Z-graph (i.e., 
pair creation from the vacuum JT3]). In what follows we adopt a Breit frame where the 
four-vector q = (q°, q) is given by q° = and q = (q x , q y , q z ) = (Q, 0, 0), while the unit vector 
n = (0,0, 1) defines the spin-quantization axis (as in Section 2). In case of the nucleon one 
has 



1*%) = F?{Q 2 ) 6, 



(22) 



leading to 



F 1 N (Q 2 



F 2 N (Q 2 



2Tr{X + } , 
2M1 



(23) 



Using Eqs. (|8|-|T3D for the nucleon wave function, the general structure of the Dirac and 
Pauli nucleon form factors (|23f) , derived from the matrix elements of the plus component of 
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the one-body e.m. current fl2~l~|), is given explicitly by: 



r 2M 

o a ,i g~ a ' 2 

(3=1,2 



E\E 2 E 3 M' q 



E\E' 2 E' 3 Mq 

^Sf wUk'^w^k^+TZ^ w* 10 (k>,p>)w 10 (k,p)+ 



+ 



fl[ n p{Q 2 ) [k\7> wi 1 (k',p')w 11 (k,p) + n ( r> w* 01 (k',p')w 01 (k,p)+ 

wl 1 {k\p')w Ql {k,p^)+T^f ) w* m (k',p')w n (k,p)}} 

where a, /? = 1, 2 and 



(a/3) 



w o(k,p) = 


1 

V2 


w s {k,p) +w s > s {k,p) 


w 01 (k,p) = 


1 


ws> a {k,P) +w A (k,p) 


w 10 (k,p) = 


1 


ws' a $,p) - w A (k,p) 


w n {k,p) = 


1 

71 


w s {k,p) - w S ' s (k,p) 



with 



dk 
dk' 



dix d& ^3^1 + 6 + 6-1] , 
dh± dk 2 ± dk 3 ± 5[ki± + k 2 _L + h±] , 

dk' 1± dk' 21 _ dk' 3± S[k' 1± - k 1± + ^q ± ] S[k' 2 ± - k 2 ± + &q±\ 
5[k' 3 ± - hx + (6 - 1) <fj 



(24) 



(25) 



(26) 



The coefficients 7^sj^s 12 , appearing in Eq. (|24]) and containing the effects of the Melosh 

rotations, are explicitly given in the Appendix B, while the quantities f^(Q 2 ) and /f (n J(Q 2 ) 
are appropriate isospin combinations of the CQ form factors, viz. 
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ftLiQ 2 ) 



1Q_ 

3 2m 



e u{d ) fr\Q 2 ) + 2e d{u) fT\Q'< 

K u{d) ff d \Q 2 )+2 Kd(u) ff u \Q 2 ) 



(27) 

wa = 0) it is 



In case of a S{7(6)-symmetric nucleon wave function (i.e., ws' a = ws> a 
straightforward to check that Eq. (El) reduces to Eqs. (10-11) of Ref. ||. 

We show now the structure of the nucleon Sachs form factors in the NR limit and in the 
ZB approximation of Ref. ||. For sake of simplicity we consider the case wa = 0, because, 
as already observed, the completely antisymmetric partial wave is very weakly coupled with 
the dominant symmetric S"-wave. Moreover, as in Ref. [§] we will limit ourselves to the case 
of flavour symmetric CQ form factors, i.e. 



f((Q 2 ) 

A A 



f(Q 2 ) 
KQ 2 ) 



(28) 



with f(Q 2 = 0) = f(Q 2 = 0) = 1. In the NR limit the Fourier transforms of the non- 
relativistic charge and magnetization densities yield (cf. Ref. ||) 



G p E {Q 2 ) 
GIAQ 2 ) 

G n E {Q 2 ) 
G n M {Q 2 ) 



F S NR (Q 2 ) + ^ R (Q 2 ) + F a NR (Q 2 ) , 

3 {Frm + ^ r (q 2 ) + « [Frm + ^ r (q 2 )} } 

F a NR (Q 2 ) + KF a NR (Q 2 )} , 



1 



-2 \ F^iQ 2 ) + ^ R (Q 2 ) + k 



F NR {Q 2, 



\^ R {Q 2 ) 



+ 



F a NR (Q 2 ) + kF?«(Q 



NRm.2\ 



(29) 



where binding effects in the nucleon mass have been neglected (i.e. M = 3m) and 



F S NR (Q 2 = \q\' z ) = f(Q 2 ) Idkdp 



Af R (Q 2 = |<?f 



f(Q 2 ) / dkdp 



w* s {kj+-q) w s (k,p) + (S -> S' s ) 
2 

w s( k > P + ^d) w s > (k, p) + {S «-> S' s ) 



F a NR (Q 2 = 151 s 



/(Q ) / dkdpw* 3 , a (k,p+-q) w s > a (k,p) 



(30) 



while F S NR (Q 2 ), A^ R (Q 2 ) and F^ R (Q 2 ) are given by Eq. (0) with f(Q 2 ) replaced by /(Q 2 ). 
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The ZB approximation can be worked out from Eq. ( Bjj) following the procedure ex- 
plained in Ref. || (cf. also Appendix B), which includes again the neglect of binding effects 
in the nucleon mass (i.e. M = 3m) for compatibility with the NR reduction. One obtains 



GUQ 2 ) 
G p M (Q 2 ) 



2M 2 ' 



3iFf(Q 2 )+Af(Q^)+K(l + 



0^- 

AM 2 ' 



f s zb (Q 2 )+^s B (Q 



F a ZB (Q 2 ) + <i + 



AM 2 



)F ZB (Q^ 



G n E (Q 2 ) 

G n M {Q 2 ) 



2M 2 



)A ZB {Q 2 ^ 



-2\f s zb (Q 2 ) + Ia zb (Q 2 ) + .(1 + ^- 2 : 



Fr{Q 2 ) + \M B {Q 2 ) 



+ 



F ZB (Q 2 ) + «(1 + 



AM 2 



)F ZB {Q 2 ) 



(31) 



where 

f zb (Q 2 ) 

A zB m 

F ZB (Q 2 ) 



f(Q 2 



1 + Q 2 /M 2 

f{Q 2 ) 
1 + Q 2 /M 2 

f(Q 2 ) 
1 + Q 2 /M 2 



[d£][dk ± ][d&±] J [w* s (k',fi) w s (k,p) + (S^ S' s ) 
[d£][dk ± ][dk' ± ] J [w* s (k',p') w s > s (k,P) + (S^S' S 
[d$[dk ± }[dk' ± } Jw* s ,(k',p') w s <(k,p) 



(32) 



with J = ^E 1 E 2 E 3 M' /E' 1 E / 2 E , 3 M . The quantities F ZB (Q 2 ), A ZB (Q 2 ) and F ZB (Q 2 ) are 

given by Eq. with f(Q 2 ) replaced by f(Q 2 ). 

From Eqs. (|2l3| - |3~2"|) it can be seen that both in the NR limit and in the ZB approximation 
G n E {Q 2 ) ^ and G P M (Q 2 ) / G n M (Q 2 ) ^ 3/2 can obtained only when the mixed-symmetry S'- 
wave is considered, more precisely through the interference between the dominant symmetric 
S'-wave and the mixed-symmetry S^-wave (cf. the definitions of A^ R (Q 2 ) and A ZB (Q 2 )). 
Such interference is expected to take place only at high values of the CQ momentum (see 
Fig. 1(b)). Finally, we stress again that in the full LF calculations (j24j) both G E (Q 2 ) ^ 
and G P M {Q 2 ) / G r l I {Q 2 ) ^ 3/2 can be obtained simply via relativistic effects only, i.e. without 
the inclusion of the effects of the mixed-symmetry S"-wave. 
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4 Results 



We have calculated the nucleon form factors adopting the nucleoli wave function of the OGE 
quark potential model of Ref. both including and excluding the effects due to the mixed- 
symmetry S"-waveP[ In the calculations we have employed only the matrix elements of the 
plus component of the e.m. current (see Eq. (|23"D ), assuming also the case of point-like 
CQ's in Eq. ( pT|) [i.e. putting f(Q 2 ) = f{Q 2 ) = 1 and k = in the right-hand side of Eq. 
(p8D l. The CQ mass m has been always taken at the value m = 0.220 GeV from the OGE 
model. Moreover, as already pointed out in Section 3, both in the NR limit and in the ZB 
approximation the binding effects in the nucleon mass have been neglected (i.e. M = 3m), 
whereas in case of the full LF calculations the nucleon mass is taken at its physical value 
(M = 0.938 GeV). 
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Figure 2. (a) The neutron charge form factor Gg(Q 2 ) vs. Q 2 within the NR limit (see Eq. (f29|)), 
assuming point-like CQ's. The data are from Ref. p5| , where the Reid Soft Core (open dots), 
Paris (full dots) and Nijemegen (diamonds) nucleon-nucleon interaction were adopted, (b) The 
ratio G P M (Q 2 )/ G^ f (Q 2 ) vs. Q 2 . The shaded area corresponds to a ±15% deviation || from the 
dipole-fit expectation G P M (Q 2 ) / G% I {Q 2 ) ~ [i v l\i n — —1.46 (dot-dashed line). Solid and dashed 
lines are the results obtained with and w/o the mixed-symmetry S"-wave, respectively. 



Our results obtained for G%(Q 2 ) and the ratio G P M (Q 2 ) / G^Q 2 ) in case of the NR limit 
(p9|), the ZB approximation fl3~T|) and the full LF calculations (0), are reported in Figs. 
2-4, respectively, and compared with the experimental data. In this respect, let us remind 

c We have obtained results very similar to those presented in this Section adopting also the nucleon wave 
function of the chiral quark potential model of Ref. ||, where the CQ mass has been chosen at the value 
m = 0.340 GeV. 
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that existing data on G P M (Q 2 ) and G^^Q 2 ) exhibit a well-known dipole behaviour, leading to 
G P M (Q 2 ) / G n M (Q 2 ) ~ Hp/fa ~ -1.46 with only a 10-=-15% uncertainty up to Q 2 ~ 1 (GeV/c) 2 
(cf, e.g., Ref. [0). 
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Figure 3. The same as in Fig. 2, but within the ZB approximation (see Eq. (|3l|)). 



From Figs. 2 and 3 it can be seen that both in the NR limit and in the ZB approximation 
the impact of the S"-wave is quite limited. As for G%(Q 2 ) (see Figs. 2(a) and 3(a)), both 
the NR and the ZB predictions are quite far from the experimental data. In particular the 
calculated neutron charge radius is still given by its canonical expectation (r)^ ~ —0.01 fm 2 
given in Table 1, which is about one order of magnitude less than the experimental value, 
( r )n,ex P = —0.113 ± 0.005 fm 2 f|. Such an underestimate is mainly due to the smallness of 
the CQ mass radii (|l9|) predicted by the OGE (as well as by the chiral) model (cf. Table 
1). Both in the NR limit and in the ZB approximation the calculated G^(Q 2 ) receives a 

El, 



positive contribution from the inclusion of the S"-wave. As pointed out in Ref. 



the 



sign of this contribution can be expected from the changes that the S"-wave produces on the 



values of CQ mass radii (|19l). As a matter of fact, the spin-spin force makes the ci-quark 
closer to the center-of-mass of the proton than the w-quark, so that (r)^ < (r)„ (cf. Table 
1). Thus the neutron charge radius receives a negative contribution (see Eq. (^3)), which 
implies an effect on G^(Q 2 ) with a positive sign. As for G P M {Q 2 ) / G^iQ 2 ) (see Figs. 2(b) 
and 3(b)), both the NR and the ZB predictions are only slightly modified by the inclusion 
of the mixed-symmetry S"-wave and therefore they are still consistent with the experimental 
data, as already noted in Ref. [Q. 

In case of the full LF calculations the effects of the S"-wave on G^(Q 2 ) are more sub- 
stantial (see Fig. 4(a)). Indeed, the LF prediction receives an important contribution from 
the inclusion of the S"-wave. The sign of this contribution is again positive as in the cases 
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of the NR limit and the ZB approximation, but with a significantly larger impact (cf. Figs. 
2(a)-4(a)). Therefore, by including the S' wave the LF result becomes much closer to the 
data. However, it should be mentioned that only ~ 65% of the experimental value of the 
neutron charge radius can be at most reproduced by our LF calculations, which, we remind, 
are based on the assumption of point-like CQ's. Note that a non-vanishing CQ size can 
affect the calculated neutron charge radius only if it is flavour dependent, i.e. different for 
u and d constituent quarks. The introduction of the effects of a possible non-vanishing CQ 
size is suggested also by the overestimate of the experimental points for Q 2 > 0.5 (GeV/c) 2 
(see Fig. 4(a)), but the estimate of such effects is beyond the aims of the present work. 




As for G P M {Q 2 ) / G^iQ 2 ) (see Fig. 4(b)), the impact of the mixed-symmetry 5"-wave is 
quite limited and therefore, at variance with the NR limit and the ZB approximation, the 
full LF calculation is not consistent with the data on the magnetic form factor ratio. In the 
next Section we illustrate that such a failure results from spurious effects due to the loss of the 
rotational covariance in the LF formalism, which occurs for approximate current operators, 
like the one-body e.m. current given by Eq. (PT|). On the contrary, we illustrate also that 
the calculations of the nucleon charge form factor G^(Q 2 ) based on the plus component of 
the e.m. current are free from unwanted spurious effects. To this end we will make use of 
the covariant LF formalism, which has been recently reviewed in Ref. [171 . 
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5 The rotational covariance problem 



As it is well known (cf., e.g., Ref. [PH), in the LF formalism the requirement of the full 
Poincare covariance of the e.m. current operator is not fulfilled by the one-body current 



This failure is related to the fact that the transverse rotations (with respect to the 
direction of the spin-quantization axis h) cannot be kinematical and therefore depend upon 
the interaction. 

An explicit manifestation of the loss of the rotational covariance is the so-called angular 
condition. As already pointed out in Section 3, the physical form factors appearing in the 
covariant decomposition of a conserved current can be expressed in terms of the matrix 
elements of only one component of the current, namely the plus component. It may occur 
however that the number of physical form factors is less than the number of the independent 
matrix elements of the plus component, obtained from the application of general symmetry 
properties to the current operator. This means that in such situations a relation among 
the matrix elements (the so-called angular condition) should occur in order to constrain 
further their number. Within the LF constituent quark model we have investigated two 
particular cases in Refs. [ 18[] and [l9j , namely the elastic form factors for the p- meson and the 



N — A(1232) transition form factors. In both cases the angular condition can be formulated 
and we have shown that the use of the one-body current (pi] ) lead to important violations 
of the angular condition, which can even totally forbid the extraction of the physical form 
factors from the matrix elements of the plus component of the current. This problem turns 
out to be particularly severe in case of "small" form factors, like e.g. the E2/M1 ratio for 



the N - A(1232) transition p 



In case of the nucleon elastic form factors it is not possible to formulate an angular 
condition, because we have the same number of physical form factors [i.e., Ff(Q 2 ) and 
F*2 '(Q 2 ) appearing in the covariant decomposition (§)] and of independent matrix elements 
of the plus component of the current (see Eq. J22])). However, this does not mean that 
rotational covariance is fulfilled. Indeed, we now make manifest the loss of the rotational 
covariance in our LF calculations by noting that the nucleon form factors can be extracted 
using not only the plus component of the e.m. current operator, but also through other 
components. As a matter of fact, adopting the Breit frame specified in Section 3, it is 
straightforward to get 

iy \n N ) = + ^p+M'iWn) (33) 

Therefore, the nucleon magnetic form factor G^(Q 2 ) can be obtained from the matrix ele- 
ments of the y-component of the e.m. current, where y is the transverse axis orthogonal to 
q±, viz. 

G N M {Q 2 ) = -^Tr[ma z ] (34) 



For the exact e.m. current the use of Eq. ( p3|) and Eq. fl34|) should yield the same result 



for Gm(Q 2 ). However, this would not be the case for an approximate current, like the one- 
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body current fl21|) employed in our calculations. As a matter of fact, the results obtained 
for G^iQ 2 ) using Eq. (p3|) and Eq. (|4]) are reported in Fig. 5, where it can clearly be 
seen that different components of the one-body current lead to quite different results. This 
finding is a clear manifestation of the presence of spurious, unphysical effects related to the 
loss of the rotational covariance in the LF formalism. It is interesting to note that we have 
carried out the calculations of G^{Q 2 ) using the y-component of the one-body e.m. current 
also within the NR limit and the ZB approximation. In both approximations the use of 
the plus and the y components yield the same result, i.e. spurious effects are found to be 
completely absent. While this finding is quite obvious in the NR limit (because the loss of 
the rotational covariance can occur only in relativistic approaches), the absence of spurious 
effects in the ZB approximation implies that the loss of rotational covariance in the LF 
formalism is related to the occurrence of the CQ initial transverse motion. Note finally from 
Fig. 5 that the LF calculations closer to the results of the NR limit are those based on the 
y-component. 




i i i I i i i I i i i I i i i 

0.2 0.4 0.6 0.8 

Q 2 (GeV/c) 2 

Figure 5. The nucleon mag netic form factor G%(Q 2 ) vs. Q 2 . Dotted lines: NR limit (Eq. (f2|)). 
Dashed lines: results obtained from the plus component of the one-body e.m. current (Eq. (p3|)). 
Solid lines: results obtained using the y- component of the one-body e.m. current (Eq. fl34|)). Upper 
and lower lines correspond to proton and neutron, respectively. Point-like CQ's are assumed. 



A possible way to circumvent the above mentioned spurious effects has been developed 
in Ref. |[20|| . There, it has been argued that the loss of the rotational covariance for an 
approximate current implies the dependence of its matrix elements upon the choice of the 
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null four- vector uj = (l,n) defining the spin-quantization axis. The particular direction 
defined by uj is irrelevant only if the current satisfies rotational covariance. Therefore, the 
covariant decomposition of an approximate current can still be done provided the four-vector 
to is explicitly considered in the construction of the relevant covariant structures. In case of 



on-shell nucleons (i.e., p 



p 



M 2 ), instead of Eq. (^|) one has 



<N(p',s')\J»(0)\N(p,s)> = u(p',s')\Ff(Q 2 )^ + F 2 N (Q 



Bi(Q) 



7 • uj 



1 



U-P M(l+T}) 



2M 



P» + B2(Q 



UJ 



uj-P 



u(p, s) 



(35) 



where F± (Q 2 ) and F 2 N (Q 2 ) are the physical form factors related to a generic approximate 
e.m. current operator J, whereas the form factors Bf(Q 2 ) (with i = 1,2,3) are unphysical 
and multiply the covariant structures depending on uj. Eventually, in Eq. (|35| ) P = (p+p')/2 

i) in the decomposition (|35|) all possible gauge-dependent 



and rj = Q 2 /AM 2 . Note that |20 



terms are forbidden because of time-reversal invariance; ii) both the physical and unphysical 
form factors do not depend explicitly upon u, because we have chosen a reference frame where 
q + = 0, and iii) for the exact current one must have Bf(Q 2 ) = 0. 
Considering the plus component in Eq. (|35D one gets 



F 1 N (Q 2 ) + ^-B^Q 2 ' 



l + rj 



Ov' N v N 



'2M 



1+71 



Bi(Q 



(36) 



which implies that for an approximate e.m. current an equation similar to (]22|) holds, but 
with the physical form factors F^(Q 2 ) and F 2 N (Q 2 ) replaced by the following combinations 



Fi{Q 2 ) 
F 2 N {Q 2 ) 



F l N (Q 2 ) + 
F 2 N (Q 2 ) + 



1 + T] 

1 



T] 



(37) 



containing the spurious form factor B±(Q 2 ). In terms of the Sachs form factors ([I]) one 
obtains 

G N E {Q 2 ) = FfiQ 2 ) - V F 2 N (Q 2 ) = F^(Q 2 ) - V F 2 N (Q 2 ) = G N E {Q 2 ) , 

G N M {Q 2 ) = F 1 N (Q 2 )+F 2 N (Q 2 )=F 1 N (Q 2 ) + F 2 N (Q 2 )+B?(Q 2 )^GZ(Q 2 ) (38) 
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Therefore, the nucleoli charge form factor G^(Q 2 ) can be safely determined from the matrix 
elements of the plus component of the e.m. current, while the same is not true for the nucleon 
magnetic form factor Gfi(Q 2 )- However, from Eq. (R3) it follows that 



jy \y»)= F?(Q 2 ) + Fi\Q 



1N1 



iQ 
2P+ 



(39) 



and therefore G^(Q 2 ) can be safely determined from the matrix elements of the ^/-component 
of the e.m. current (as in Eq. ([34])). 

In Fig. 6 we have reported our final results based on the use of Eq. ( [23] ) for G%{Q 2 ) 
and on Eq. (]34] ) for the ratio G P M {Q 2 ) / G^iQ 2 ) . It can be seen that, provided spurious 
effects are properly avoided, our LF results are fully consistent with the experimental data 
on G p M (Q 2 )/GUQ 2 )- 
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Figure 6. The same as in Fig. 2, but within the LF approach using the plus component of 
the one-body e.m. current (see Eq. (|23|)) for the determination of G 1 ^{Q 2 ) and the transverse 
y-component (see Eq. (p3)) for the calculation of the nucleon magnetic form factors. 



To sum up, the most relevant SU(6) breaking exhibited by the experimental data, namely 
Ge(Q 2 ) 7^ and G P M (Q 2 ) / 'G\ { (Q 2 ) ^ —3/2, can be understood qualitatively as well as 
quantitatively within the CQ model, provided the effects of the spin- dependent components 
of the effective quark-quark interaction are taken into account as well as the relativistic 
effects arising from the LF composition of the CQ spins are properly considered. 
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6 The ratio G P E (Q 2 )/G P M (Q 2 ) 

In this Section we will briefly address the issue of the interpretation of the recently observed 
PH deviation of the ratio ix p G p E {Q 2 ) / G P M {Q 2 ) from the dipole-fit expectation fi p G p E {Q 2 )/ 
G P M {Q 2 ) = 1, where n v = G P M (0) is the proton magnetic moment. First of all let us note that 
if the mixed-symmetry S"-wave is neglected and point-like CQ's are assumed, the NR limit 
( p9|) predicts a ratio n P G p E {Q 2 ) / G P M {Q 2 ) = 1 (i.e., coinciding with the dipole-fit expectation), 
while the ZB approximation (||) yields fj, p G p E (Q 2 )/G p M (Q 2 ) = 1 - Q 2 /2M 2 , where now \i v 
is the calculated proton magnetic moment. 
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Figure 7. Ratio ^ P G P E (Q 2 )/G P M (Q 2 
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limit ( 

(IT 



Q 2 (GeV/c) 2 



versus Q 2 . Open squares are the recent experimental points 
[10 1- ( a ) Dotted, dashed and solid lines correspond to the predictions of the NR 
and the ZB approximation (pip, as well as to our full LF calculations [namely Eq. 
for G P E (Q 2 ) and Eq. (|34|) for G P M (Q 2 )], respectively. Thick and thin lines are the results 
obtained with and without the effects from the mixed-symmetry S"-wave (dotted and dashed thin 
lines are undistinguishable from the thick ones), respectively, (b) The thick solid line is the result 
obtained using the full OGE interaction model @ (yielding (p±) = 0.58 GeV), while the dashed 
line corresponds to the case of its linear confinement term only (corresponding to (p±) = 0.33 GeV). 
The thin solid line corresponds to the use of the chiral quark potential model of Ref. [||] (having 
(p±) = 0.61 GeV), while the dotted line is the result obtained adopting a pure S'C/(6)-symmetric 
HO ansatz corresponding to (p±) = 1.0 GeV, respectively. In all the calculations point-like CQ's 
are assumed. 



The full LF calculations, based on Eq. (0) for G P E (Q 2 ) and Eq. fl3|) for G P M (Q 2 ), and 
the results obtained in the NR limit and the ZB approximation are reported in Fig. 7(a) 
and compared with the recent J LAB data ||10|| . It can be seen that: i) the impact of the 
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mixed-symmetry S"-wave is very limited in each approach, and ii) the predictions of both 
the NR limit and the ZB approximation are completely at variance with the JLAB data; 
in particular, the ZB approximation predicts a negative value for fi p G E (Q 2 )/G p VI (Q 2 ) for 
Q 2 > 2M 2 (with M = 3 • 0.220 GeV for the NR and ZB calculations); this clearly signals 
that the applicability of the ZB approximation should be limited only to low values of Q 2 , 
namely Q 2 « M 2 (cf. also Ref. ||). To sum up, our findings imply that a suppression 
of the ratio ii p G p E {Q 2 ) / G P M {Q 2 ) from unity can be expected in the CQ model provided the 
relativistic effects generated by the Melosh rotations of the CQ spins are taken into account. 

The LF results obtained adopting the OGE and chiral quark potential models are re- 
ported in Fig. 7(b) and compared with calculations corresponding to various values of the 
average CQ transverse momentum (p±). It can clearly be seen that the calculated suppres- 
sion of the ratio /i p G E (Q 2 ) /G P M (Q 2 ) from unity exhibits a slight dependence on the value of 
(p±), and moreover the results corresponding to the OGE and chiral quark potential models 
are very similar and compare quite favourably against the recent JLAB data. A nice re- 
production of the latter appears to be achieved when (p±) ~ 1 GeV [see dotted line in Fig. 
7(b)], but it should be reminded that our results have been obtained assuming point-like 
CQ's, while the introduction of CQ form factors is expected to affect (at least partially) the 
ratio ^ P G p E {Q 2 )/G p M {Q 2 ). 



7 Conclusions 

In this contribution we have investigated the effects of both kinematical and dynamical SU(6) 
breaking on the nucleon elastic electromagnetic form factors within the constituent quark 
model formulated on the light-front. We have shown that the most relevant SU(6) breaking 
exhibited by the experimental data, namely G E {Q 2 ) ^ and G P M (Q 2 ) / G^iQ 2 ) ^ —3/2, can 
be understood qualitatively as well as quantitatively within the constituent quark model, 
provided the effects of the spin-dependent components of the quark-quark interaction are 
taken into account as well as the relativistic effects arising from the light-front composition 
of the constituent spins are properly considered. We have indeed shown that the evaluation 
of the nucleon magnetic form factors has to be performed using the transverse ^-component 
of the e.m. current in order to avoid spurious, unphysical effects related to the loss of the 
rotational covariance in the light-front formalism, while the nucleon charge form factor can 
be safely extracted from the matrix elements of the plus component of the current. 

Finally, we have shown that a suppression of the ratio G P E (Q 2 ) / G P M {Q 2 ) with respect 
to the dipole-fit prediction can be expected in the constituent quark model provided the 
relativistic effects generated by the Melosh rotations of the constituent spins are taken into 
account. The strength of the suppression exhibits a slight dependence on the value of the 
average quark transverse momentum in the nucleon; moreover, the results of the calcula- 
tions based on the nucleon wave functions arising from two of the most sophisticated quark 
potential models |7|, || , compare quite favourably against the recent data (10| from Jefferson 
Lab. 
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Appendix 

A. Expansion of the wave function in the harmonic oscillator basis 
having definite per mutational symmetry 

The wave functions ws(k,p), Ws> 3 (k,p), ws> a (k,p) and WA(k,p), appearing in Eq. (|TT|), can 
be constructed by expanding them onto the complete harmonic oscillator {HO) basis, viz. 



w [f] {k,p) 



lh (HO,[f]) 



(k,p) 



(40) 



(HO \ f]) ~ * 

where Wp^M L \k,p) is a function of the HO basis corresponding to total orbital angular 
momentum L and its projection Ml, p stands for all the other quantum numbers necessary 
to identify the basis functions (see later on) and [/] identifies the component of the nucleon 
wave function (i.e., [/] = S, S' s , S' a , A). Finally, the quantities are variational coefficients, 
that can be determined by applying the Raleigh-Ritz variational principle to the Hamiltonian 
of the three-quark system. 

The basis function LJ ^ L (k,p) can be written in the form 



(HO, {f]) 
p,LM L 



(HO,) 



(41) 



where p is a shorthand notation for the HO radial and orbital quantum numbers p 
{nk,£k,n p ,£ p } and ^pLMr i^iP) ls the usua l HO wave function given explicitly by 



l{HO) 
J p,LM L 



'P,L,M h 

(k,p) = R { n H k ° ] (k) R { ™\p) (ikm k J pi mJLM L )Y ekmk (k)Y epmp (p) (42) 



m k m v 



with B}^ k f k \k) and R^f\p) being the radial HO functions depending on the HO length 
duo [namely, taking into account the definitions (|T2]), one has: R^ ?Mk) oc exp(—k 2 /2a 2 HO ) 

and R^f\p) cx exp(—3p 2 /8a 2 HO )}. Thanks to the use of only one HO length, in Eq. ( |4l"l) 
the sum over p is limited to the HO functions having the same number of HO excitation 
quanta Nho = N k + N p = In^k + 2n p + £ p , the same orbital parity n = (— ) £k+e p and the 
same orbital angular momentum L and its projection Mi. In Eq. ([|1|) the matrix U ~ , 
which is independent on M^ thanks to the Wigner-Eckart theorem, is explicitly given by 



U. 



L,[f]=S 



PP 



1 + 



\ 1 + 5h k n p 3(> k i p 



ep (p;L\p; L)A dfi 



U 



L,[f]=S' a 
PP 



u. 



L,[f]=S' a 
PP 



V2 



V2 



1 + 



1 " R' 



u. 



L,[f]=A 



PP 



^+^(p-L\p;L)[A dA + A d , 2 ] , 
(p;L\p; L) [-A djl + A d>2 ] 
(p;L\p; L)A dfi 



(43) 
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where (p; L\p; L) is a shorthand notation for the Brody-Moshinsky coefficients [^T], d 
N k -N p = 2n k + 1 



2n r 



and A d)m = 1 if d = m [mod(3)] while A d m = otherwise. 



{n fc ,£ 



}, but the following constraints apply: 



Finally, the index p stands for p 
Uk > n p , ii) £k > £ p if rik = n P , and iii) hk ^ n p for [/] = A and [/] = S if L is odd. 

Note that in the right-hand sides of Eq. fl4*3| ) the presence of the factor (— ) ip is related 
to the definition of the jacobian variable p as given in Eq. (|i~2l). 

In the calculations presented in this work we have employed the HO basis up to 20 
HO excitation quanta (i.e. Nho < 20), which corresponds to a total of 67 basis states 
for the symmetric S'-wave, 94 states for the mixed-symmetry S"-wave and 31 states for 
the antisymmetric A-w&ve. By combining these configuration states with the proper spin- 
isospin functions, one gets a total of 192 (completely symmetric) HO basis states for the 
expansion of the (canonical) nucleon wave function ( 11 
functions 

symmetry, is 572 for Nho — 20 and L = Ml = 



Note that the number of HO 
p H LAi L (k,p), which appear in Eq. (^TJ) and do not posses any definite permutational 



B. Evaluation of the coefficients a 



In Eq. (|24D the coefficients Ttp s contains the effects of the Melosh rotations of the CQ 
spins and are defined as 



S' 12 S 12 



n { T 2 /3) 

S' 12 Sl2 



V N 

i 



, E ((^2 I) \vn\ ^ 0? U | (s l2 ±) \v' N ) (\v' N \ a y \\v N ) (44) 

vnv' n 



where V) is given by Eqs. (|9|-|T0D, O 1 is the identity 2x2 matrix and O 2 
the spin of the particle 3, and 

S12 ~) g^Jv) = ^(~vi~V2\Si2M s ) (Si 2 M s -v 3 \-v N ) 

M s 



-ia y acting on 
(45) 



The sums over the spin projections yield 



/c oo 



K 



(12) 
00 



-r>(21) 
/c 00 



(22) 
00 



Af 
M 
J\f 
M 



A X A 2 + B l ■ B 2 \ A 3 , 
A X A 2 + B x ■ B 2 ] A 3 , 
A X A 2 + Bt ■ B 2 ] B 3 , 
A X A 2 + B 1 ■ B 2 ] B 3 , 



(46) 
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(ii) 

01 



-r>(12) 

/Vfi 



-01 



n, 



(21) 
01 



-,(22) 
'^01 



A 2j B! - AiB 2 - fix x 5 2 j • B 3 , 
A 2 B 1 - A X B 2 -B x x B 2 ] ■ B 3 , 



-^A/" — A ± B 2y + B 2z Bi x — B 2x B lz ) A 3 + 

(AiB 2x — A 2 B lx — B 2y B\ z + B 2z Bi y ) B 3z — 
{A\B 2z — A 2 B\ Z — B 2x B\ y + B 2y B\ x ) B 3x ] , 

— j=Af (A 2 Bi y — A\B 2y + B 2z B lx — B 2x B lz ) A 3 + 



{A x B 2x — A 2 B lx — B 2y B lz + B 2z B ly ) B 3z — 
{A\B 2z — A 2 B\ Z — B 2x Bi y + B 2y Bi x ) B 3x 



(47) 



K 



(ii) 

10 



-77(12) 
'Mo 



K 



(21) 
10 



-77(22) 
'MO 



j=N \A 2 B X - A X B 2 + 5 lX B 2 ] ■ B 3 , 

V3 L J 



— j=N 



A 2 B X - A X B 2 + B l xB 2 



B, 



v/3 



M [(A 2 Bi y — A\B 2y — B 2z B\ x + B 2x Bi z ) A 3 - 



{A\B 2x — A 2 B\ X + B 2y B lz — B 2z B ly ) B 3z + 
{A\B 2z — A 2 Bi z + B 2x B ly — B 2y B lx ) B 3x ] , 
1 



M (A 2 B ly — AxB 2y — B 2z B lx + B 2x B lz ) A 3 - 



{A\B 2x — A 2 Bi x + B 2y B lz — B 2z B ly ) B 3z + 
{A\B 2z — A 2 B lz + B 2x B ly — B 2y B lx ) B 3x 

= \m[(?>A 1 A 2 -B 1 -B 2 )a 3 + 2(a 1 B 2 + A 2 B 1 )-B 3 



-77(12) 

'Ml 



-M 



(SA 1 A 2 - B x ■ B 2 ) A 3 + 2 (A X B 2 + A 2 B^ ■ B 3 



(48) 
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-p (21) 
'Ml 



-77(22) 
'Ml 



2 [B lv B 2 + B 2v B! ■ 5 3 - (A 1 A 2 + B x ■ B 2 ) B 3y + 



2 {A 1 B 2y + A 2 B ly ) A 3 ] 
1 



2 B lv B 2 + 5 2v 5i ■ 5 3 - A X A 2 + B ± -B 2 ) B 3y + 



where 



Af 



A, 



-J\f 



2 (A 1 B 2y + A 2 B ly ) A z 



3 

n 



(49) 



\J(m + ^M o y + kl 



m + i 3 M' Q f + k>) 



j-L 



B 



B jz 

B 3x 

B 3 y 

B 3z 



(m + ijM'o) (m + ^Mq) + P i± ■ k j± , 
(m + &M' )k 3x - (m + £ 3 M )£;' 3a; , 
(m + ^M )k' jy - (m + CjM' )k jy , 
(to + ijM' )k jx - (to + )//,.-■ , 



. t . — h' ■ k- 

J x 3V 3V 3 X ' 



& 3i&3i/ + k' 3y k 3x , 

-(m + £ 3 M' )(m + £ 3 M ) - /c'^/c^ + fc^A^ 
(to + £ 3 M' )k 3y + (to + £ 3 M )£;' 32/ 



(50) 



The approximation of Ref. [0] (see Eq. ([Jl])) can be obtained from the full 
results ( p4] , f46| - [4^ ) making the assumptions explained in Ref. ||. One easily obtains: 1Zq{^ = 



K 



(a/3) 
10 



and 



/c 00 



m 



K 



(12) 
00 



m 2 + Q 2 /9 
1 Q 



m 



11, 



(21) 
00 



3 2m y/ m 2 + Q2/9 
2 Q 



m 



3 2m / m 2 + Q2/9 



-r>(22) 
/c 00 



m 



m 2 + Q 2 /9 



(51) 
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K 



(ii) 
11 



4m 2 - 3Q 2 /9 m 
Am? + Q 2 /9 ^ m 2 + Q2/g ' 



^(12) 1 Q m 



9 2m / m 2 + g 2 /g 



(21) 2 Q 4m 2 - Q 2 /9 m 

'Mi 



3 2m4m 2 + Q 2 /9^/ m 2 + g2 /9 ' 

= \ , m • (52) 

3 + Q 2 /9 
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